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Abstract 

According to self-similarity hypothesis, the thermodynamic limit could be defined from the scaling laws for the system self- 
similarity by using the microcanonical ensemble. This analysis for selfgravitating systems yields the following thermodynamic 
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limit: TV — » oo keeping constant E/N~s and LTV 3, in which is ensured the extensivity of the Boltzmann entropy Sb = 
In W (E, TV). It is shown how the consideration of this thermodynamic limit allows us to explain the origin of dynamical 
anomalies in numerical simulations of selfgravitating systems. 

PACS numbers: 05.20.-y; 05.70.-a 



I. INTRODUCTION 

In the last years there is a special interest to perform 
a well-defined thermodynamical description for selfgrav- 
itating systems and many interesting results have been 
obtained 0, S S S S H S H Eg . In ref.[ll|, de Vega 
and Sanchez have pointed out that a kind of thermody- 
namical limit of a self-gravitating system can be defined 
if one considers what they call the diluted limit: send the 
number of particles, TV, and the volume, V , to infinity, 
keeping constant the ratio TV/ V 1 ' 3 instead of the density, 
N/V. Recently, V. Laliena in ref.^2 showed that this 
thermodynamic limit suffers from the same problems as 
the usual thermodynamical limit and leads to divergent 
thermodynamical functions. 

In the recent paper , we proposed a new thermody- 
namic limit for selfgravitating system which is obtained 
from the self- similarity properties [lH [ThL Hi! ITtI ITsj of 
this kind of systems: 



E i 
TV — » oo, keeping constant — =- and LTV 3 . (1) 

TVs 

As already shown in ref.EH, the consideration of this 
thermodynamic limit allows us to perform a well-defined 
thermodynamical description for selfgravitating systems 
by using almost the same arguments used by the stan- 
dard Thermodynamics for describing the extensive sys- 
tems. The aim of this paper is to show that the dynami- 
cal anomalies in numerical simulations of selfgravitating 
systems disappear when the self-similarity properties of 
this kind of systems are taken into account. 
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II. WHAT IS THE SYSTEM SELF-SIMILARITY? 

Self-similarity is a natural generalization of the exten- 
sivity of the traditional system, which should be taken 
into account in order to perform a well-defined thermo- 
dynamical description based on the microcanonical en- 
semble [HE EH. Basically, it can be obtained from 
the asymptotical scaling behavior of thermodynamical 
variables and the microcanonical phase-space volume W 
when the many particle limit TV — > oo is invoked. 

There are many studies in long-range Hamiltonian sys- 
tems where the validity of the thermodynamical descrip- 
tion in the thermodynamic limit is intimately related to 
the scaling with TV of the thermodynamic variables and 
potentials. The problems with the nonextensive nature of 
those systems are avoided using the physically unjustified 
Kac prescription |l9j| in which the coupling constants are 
scaled by some power of TV in order to deal with an ex- 
tensive total energy E. Self-similarity can not be reduced 
to the Kac prescription, because the energy extensivity is 
not demanded in the self-similarity framework. Geomet- 
rical aspects of the probabilistic distribution function of 
the statistical ensembles allows us to define a well-defined 
thermodynamic formalism even in cases where the exten- 
sive nature of systems can not be ensured 0, 0, . 

The type of self-similarity scaling laws determines 
which is the generalized Boltzmann entropy [lj which 
should be considered in the application in order to guar- 
antee the ensemble equivalence of the microcanonical en- 
semble with some adequate generalization of the canon- 
ical ensemble |l4j . Ordinarily, scaling laws are exponen- 
tial, it means that the microcanonical accessible volume 
W has an exponential growing with TV in the thermody- 
namic limit. There is not a priory reason for supposing 
that a Hamiltonian system should necessarily satisfy this 
kind of growing in the thermodynamic limit, although 
there is a huge world of Hamiltonian systems that do 
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obey it. In ref.|l6( we used the self-similarity hypothesis 
in order to find the necessary conditions for the validity 
of Tsallis' Statistics |2£| by considering these arguments 
and using a procedure inspired in the one used by Gross 
in deriving his microcanonical thermostatistical theory 
|2l| . It was proved in that paper the nonuniqueness of a 
thermodynamic formalism based on the microcanonical 
ensemble. In this sense, self-similarity also difers from 
the Kac prescription. 

III. THE SELF-SIMILARITY SCALING LAWS 

In this section is recalled how can be obtained the 
thermodynamic limit by considering the system self- 
similarity properties for the following tridimensional self- 
gravitating Hamiltonian system: 



Thus, when the thermodynamic limit is taken into 
account, the Boltzmann entropy is extensive: 

S B =lnW~N. (8) 

There are some connections suggesting the validity of 
the thermodynamic limit (Q. It was proposed in ref . jig . 
an alternative model for the classic isothermal model of 
Antonov p| which uses an energetic prescription instead 
of a box renormalization in order to avoid the long-range 
singularity of the A-body selfgravitating systems. In that 
paper, a classical gas of identical particles with mass m 
was considered. Taking into account the rest energy of 
the particles, the nonrelativistic limit is valid when the 
absolute value of the mechanical energy of the system is 
much smaller than its rest energy : 



H = T + U 



N 1 

fc=i 



N 

E 

j>k=i 



<t>{rjk) 



(2) 



Here — |ij — rj.| is the distance between the j-th and 
k-th particles, being <p (r) the selfgravitating potential 
which behaves as Newtonian interaction when r — ► oo: 



<j) (r) ~ 1/r. 



(3) 



It is easy to show that the microcanonical phase-space 
accessible volume W is given by: 



e e($ ) 



< mc 2 N, 



(9) 



where £q is the characteristic energy of this model: 



e = 2ir- 



H 2 ' 

However, this condition can not be satisfied for an ar- 
bitrary number of particles. In fact, when N tends to 
N : 



W 



NT 



i r N 

T[ d h k [E 

{—) Jv N k=l 



(4) 



where V is the tridimensional volume in which the system 
is enclosed. Since cj> (r) describes a long-range interaction, 
the potential energy U for N large could be estimated as 
follows: 



U ~ N 2 L- 



(5) 



where L is the characteristic linear dimension of the sys- 
tem. The factor N 2 takes into account the contribution 
of all couples of particles, while the factor L _1 appears 
as consequence of the behavior of the potential <j> (r) Q 
for r large . The total energy E will exhibit an identical 
behavior: 



E - N 2 L- 



(6) 



Thus, when N is sent to infinity, the microcanonical 
phase-space accessible volume W depends on N and L 
as follows: 



/ mc 2 \ 4 




I 2 '-) " 


(1) 
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which corresponds to a characteristic mass Mq: 



'fic\ 5 1 
M = N m=[-) - 2 , 



(10) 



(11) 



the model loses its validity. Everybody can recognize the 
fundamental constant of the stellar systems, which has 
much to do with the stability conditions of the stars (see 
in ref. 22])- Note that this constant appears as a con- 
sequence of assuming the energy ./V-dependence (Q, so 
that, it could not have been obtained if another thermo- 
dynamic limit had been adopted. A consequent analysis 
of these massive systems should be performed taking into 
account the relativistic effects. 

Another link is found in the well-known white dwarfs 
model based on the consideration of the Thomas-Fermi 
method to describe the state equation of the degenerate 
nonrelativistic electronic gas, whose pressure supports 
the hydrostatic equilibrium of the star: 



W 



N\T 



' 3N 



- [N 2 L- X ] ^ L 3N 



LNz 



(7) 



A(j> = -AirGp, 



(12) 
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where <fi is the Newtonian potential, being p the mass 
density: 



2 s mi 



P = ^ ^ (13) 

where p, is the number of nucleons per electron, m and 
m e are the nucleon and electron masses, being <f>s the 
potential at the star surface. From this model it can be 
easily derived the iV-dependence of total energy and lin- 
ear dimension of the system given in Q by using a sim- 
ple dimensional analysis. This coincidence is not casual: 
for the nonrelativistic particles the thermodynamic limit 
only depends on the dimension of the physical space. 

According to the exposed above, during the many par- 
ticle limit, N — > oo, the selfgravitating Hamiltonian sys- 
tems exhibit self- similarity under the following scal- 
ing laws: 

N N (a) = aN , \ 

E -> E (a) = otl E a , \ W (1) -> W (a) = T [W (1) 
L — > L (a) = a~zL , ) 

(14) 

where the functional T [W, a] defines an exponential scal- 
ing laws: 



F[W,a] =exp [ahxW] 



(15) 



being 

W (a) = W [E (a) , N (a) , L (a)} . (16) 
As already mentioned in our previous works 

mill 

llH lri| . a system exhibits self-similarity under certain 
scaling transformations of its fundamental macroscopic 
observables when the functional T \W, a] obeys to the 
general condition of self-similarity [lj|: 



T \T [a:, a{\ ,a 2 ] = J 7 [x, a\a 2 ] 



(17) 



Only in this case, the scaling transformations constitute 
a group of transformations compatible with the system 
macroscopic properties, whose ordering information is 
contained in the microcanonical accessible phase space 
volume W . Wc suppose that when the thermodynamic 
limit is carried out by using these special scaling laws: 

• all those system self-similarity properties are pro- 
tected, even though, its dynamic characteristics, 

• as well as under certain conditions, the microcanon- 
ical ensemble could be equivalent to certain gener- 
alization of canonical ensemble. 

Contrary, we guess that the non consideration of the 
self-similarity scaling laws should provoke a trivial en- 
semble inequivalence, as well as others anomalous behav- 
iors According to our viewpoint, the Boltzmann's 
definition of entropy is applicable to this system, which 
allows us to classify the selfgravitating model as a pseu- 
doextensive system [l5| . 



IV. DISCUSSIONS AND CONCLUSIONS 

Many investigators do not pay so much attention to the 
analysis of the system self-similarity when it is performed 
its macroscopic description. However, the non or bad 
consideration of the system self-similarity scaling laws 
leads in many cases to the trivial nonequivalence of the 
statistics ensembles as well as other anomalous behaviors, 
such as the reduction of the mixing and slow relaxation 
regime towards the equilibrium. 

As an example of the above affirmation are the anoma- 
lies presented in the dynamical study of the selfgravitat- 
ing systems performed by Cerruti-Sola & Pettini in the 
rcf. 23] . In that paper, they observed a weakening of the 
system chaotic behavior with the increasing of the parti- 
cles number N . It is very well-known the consequence of 
this fact on the ergodicity of the system: the chaotic dy- 
namics provides the mixing property in the phase space 
necessary for obtaining the equilibrium. In that exam- 
ple the chaoticity time grows with N, and therefore, the 
. oi\ systems could expend so much time to arrive to the equi- 
librium. 

In a recent paper p4|. Taruya and Sakagami carried 
out a iV-body numerical experiment of the same situa- 
tion as investigated in classic papers [5j,|6j. That is, they 
confined the N particles interacting via Newton grav- 
ity in a spherical adiabatic wall, which reverses the ra- 
dial components of the velocity if the particle reaches 
the wall. The typical timescales appearing in this sys- 
tem are the free-fall time, Tff = (Gp)^ 1 ^ 2 , and the 
global relaxation time driven by the two-body encounter, 
Treiax = (0.1 iV/ In N)T ff [2j|. Since they considered the 
density p as A-independent quantity, the relaxation time 
T re iax diverges when the thermodynamic limit is invoked. 

These facts can be interpreted as the non- 
commutativity of the thermodynamic limit (N — > oo) 
with the infinite time limit ( t —> oo): when the first is 
performed before the second, the systems will not relax 
to the microcanonical ensemble. However, our analysis 
allows us to understand the origin of all these anomalous 
behaviors. 

It was considered in the study [23| that the energy is 
scaled proportional to N during the realization of the 
thermodynamic limit, which does not take into account 
the system self-similarity. It is very easy to see that the 
instability exponent Xh versus energy dependence (Fig- 
ure 4 of the ref. 23]) is corrected when the self-similarity 
scaling law for the energy is considered: its proportional- 
ity to N~s. The weakening of chaoticity disappears when 
self-similarity is taken into account. Anyway, in spite 
of they assumed a different scaling law for the energy, 
they obtained the correct dependency of the instability 
exponent Xh with the energy per particle: the power law 

The dynamical anomalies in the Taruya and Sakagami 
analysis are also explained by the non consideration of 
self-similarity. When the thermodynamic limit Q is 
taken into consideration, system volume scales as ~ 1/A, 
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and therefore, the characteristic density p scales as ~ TV 2 . 
Thus, Tff ~ l/N, and hence, T relax ~ l/\nN -> 0. 
So that, the origin of the dynamical anomalies is a triv- 
ial consequence of the consideration of a thermodynamic 
limit incompatible with the system self-similarity prop- 
erties of the self-gravitating systems. 

There are some other examples of dynamical anoma- 



lies which at first glance could lead to the inapplicability 
of the microcanonical or the canonical ensemble (see for 
example in refs. [H El S H3, El HI HI IP ) ■ However, 
it could be proved that all these dynamical anomalies 
can be explained by the same arguments exposed in this 
paper. 
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